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Abstract 

We investigate the possibility to disentangle the radions in the Randall-Sundrum scenario 
from higgs boson at the next generation high energy e + e _ linear colliders. Due to trace anomaly, 
the radion coupling (and in turn the branching ratio) to gluons is enhanced over the same for the 
Standard Model (SM) higgs. We study the radion <f> production at electron-positron colliders, 
via the process e + e~ — > w<j), and propose to investigate radion decay to a pair of gluons. 
At 500 GeV center of mass energy, our signal is not very promising. We find our results are 
encouraging for finding radion and differentiating it from higgs for center of mass energies of 1 
TeV and above. 

Recently proposed models [1, 2], aiming to cure the hierarchy problem between electroweak (EW) 
and Planck scale, offer an exciting possibility to test the gravitational interactions at the TeV 
colliders. The models in [1, 2] require our universe to have (3 + n) + 1 space-time dimensions, of 
which the extra, n, space-like dimensions are compactified. The ADD [1] model requires relatively 
large compactification radius (~ 1 mm). In the Randall-Sundrum (RS) scenario [2], which is 
the focus of this analysis, there is only one extra space-like dimension compactified on an S1/Z2 
orbifold. Unlike the ADD, RS-scenario does not require a large compactification radius for the 
extra compactified space-like dimension. The radius of compactification is of the order of the 
Planck length, and interestingly it is a dynamical object. Moreover, the assumed geometry of 
space-time is non-factorisable in the sense that 3 + 1 dimensional metric is scaled by an exponential 
warp-factor depending on the extra space-like dimension. Apart from the Kaluza-Klein tower of 
gravitons, the low energy effective theory has also a graviscalar, called radion. Goldberger and 
Wise [3] have proposed a method to generate a potential for this scalar, by introducing a scalar 
field in full five dimensional manifold. This in turn dynamically generates a vacuum expectation 
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value (VEV), A^, for the radion. Following the proposal in [3], A^, comes out to be of the order of 
TeV without finetuning the parameters. The Kaluza-Klein excitations of the bulk fields are of the 
order of a few times TeV [3]. Assuming a stabilization similar to [3], the radion is likely to be the 
lowest lying gravitational state. Here we will assume that only radion is of interest in the studied 
experimental situations. 

Radion <p couples to the Standard Model particles in a model independent fashion via the trace 
of the energy momentum tensor, 

Lint = ~T~ TV 0. (1) 

This implies that radion interaction with SM fields is very similar to the higgs boson, but with a 
suppressed strength depending on the value of A^. Apart from eq. (1), radion can have a mixing 
with the SM higgs via the following term in the action: 

S = -tJ d A x^-R(g vis )H^H. (2) 

Here the Ricci scalar R (g v i s ) corresponds to the induced four dimensional metric, g v i s , on the 
visible brane, H is the electroweak higgs boson, and £ is the mixing parameter. 

Phenomenology of radion has been studied extensively in recent literature in the context of 
collider experiments like LHC, a linear collider, or in the cases like muon (g — 2), K° — K° mixing, 
and electroweak p-parameter [7]. Effects of radion have also been studied in the context of unitar- 
ity in gauge boson scattering [8]. While the low energy experiments try to constrain the radion 
parameters, collider studies are naturally more focussed on the possible search strategies of this 
particle at the planned experiments. From the collider studies it is evident that LHC would be the 
best place to discover /exclude this kind of a scalar if the radion vev is in the ball park of a TeV or 
so. 

At the same time, in order to make correct theoretical deductions, it is crucial to know the 
real identity of a scalar detected in an experiment. In this work we consider the detection and 
identification of radion, which may be difficult, since the decay modes of radion and higgs are 
identical. A recent study [9] aimed at the high energy e + e~ collider, proposes associated production 
of a radion with a SM higgs mediated by KK gravitons. This particular final state is an outcome of 
higgs-radion mixing. Higgs-radion mixing have also been investigated to study the complementarity 
of higgs/radion signals at e + e~ colliders [10]. In the following we will study the separation of the 
radion signal from the higgs signal at a linear collider, with varying radion-higgs mixing parameter. 

As mentioned earlier, radion couples to the SM particles via the trace of energy momentum 
tensor. Being massless, gluons and photons cannot couple to a radion at classical level. At quantum 
level they do couple to radion, with a coupling proportional to the /3-function coefficients of U(l)y 
and SU(3)c, respectively. It turns out that for A^ ~ 1 TeV, radion couplings to a pair of gluons 
or photons are enhanced over the corresponding couplings of higgs. This was the key feature of all 
the previous exercises looking for a radion at hadronic machines. At an e + e~ collider, we propose 
to produce radions in the process e + e _ — > vucj), and to study radion decay to gluon gluon. Though 



radion branching ratio to gluon gluon is quite sizeable, this decay channel cannot be used at hadron 
colliders due to overwhelmingly large QCD background. At an e + e" machine, one can hope to cope 
with two jets with missing energy-momentum final state. 

After a brief discussion of higgs-radion mixing, we will continue with the signal and background 
analysis. The action in eq. (2) leads to a term bi-linear in curvature and higgs in the Lagrangian 
[11, 12], 

C = -6£ft 2 (□ In Q + (V In ft) 2 ) H*H, 

where 

n = e -(7/«)^) > 7 = U /A . 

Clearly, this induces a kinetic mixing between higgs and the radion. Additional mixing is introduced 
by the fact that both the neutral component of the higgs (h) as well as the radion ((f)) acquire 
VEVs ((h) = v). To obtain fields with canonical quantization rules, it is necessary to make field 
redefinitions: 

, I d>\ ( cos 9 — sin 9 \ 

= Z R M- 1 M=[ (3) 

\h J \ ZfisinO — 6^7 cos 6* Zr cos 9 + 6^7 sin 9) 

= 1-6^(1 + 60 a»d tanM = - p igg^_ i (4) 

where irih and are the higgs and radion mass parameters in the Lagrangian. In the limit £ — ► 0, 
we recover back the SM higgs from h' . 

Requiring that the kinetic terms for the physical fields h! and <p' be positive, restricts us to [12]: 

^ (1 + \A + w) <^ <^(yr+w-i). (5) 

For A^ = 1 TeV, this translates to —0.75 < £ < 0.59. As seen from eq. (3), the mixing matrix of 
radion and higgs is not unitary. Therefore it is not always straightforward, which particle should 
be called higgs and which should be called radion. We will always call <fi' radion and h' higgs in 
the following calculations. 

The above redefinition of the fields leads to the following interactions between scalars and 
fermions or gauge bosons: 




^<S>u' , a! 
a34—n + oi2<?> 



(6) 



where 034 and 012 are defined in terms of the elements of the matrix M: au = (M-n + M.21/1) 
and a 34 = ^ (M22 + 7X12). 

The effects of mixing on the higgs and radion couplings to fermions and gauge bosons are 
presented in fig. 1, see eq. (6) for the definitions of (034/1;) and (a^/A^,). From the figure it is 
evident that apart from the region around the special value of £ = |, radion couplings to gauge 
bosons or fermions are bigger than the respective higgs couplings. Far enough from |£| = |, cp' 
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Figure 1: Comparison of radion (marked <f>') and higgs (marked h') couplings to the fermions or 
gauge bosons as denned in eq. (6). The used mass parameters are rrih = 150 GeV and = 250, 500 
GeV. 

couplings are enhanced over the h! couplings. This phenomenon is purely due to the nature of the 
higgs-radion mixing. While calculating the radion (higgs) couplings, we have fixed = 150 GeV 
and = 250 or 500 GeV. 

The couplings to gluons can be written as follows [13]: 
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where 63 is the QCD /3-function coefficient. i^/ 2 is the form factor from (heavy quark) loop effects. 
In each of these couplings the first term proportional to 63 is coming from the trace anomaly. We can 
see from eq. (7) that the vertices higgs/radion -gluon-gluon have new contributions, which change 
the production and decay of the higgs boson. Couplings of radion/higgs to a pair of photons can 
also have anomalous contributions which we do not write here explicitly and which can be found 
elsewhere [13]. We do not write either the trilinear couplings involving radion and a pair of higgs 
or vice-versa, as we are not interested in these couplings in this work. 

Without going into further details of the gluon couplings, we choose to plot in fig. 2 the radion 
branching ratio to gluons for both vanishing and non-zero values of £. For purpose of comparison 
we also present the same branching ratios for the higgs boson. For £ = 0, — > gg branching ratio 
is almost two orders of magnitude higher than the same for h. The sharp fall of the branching 
ratio around 160 GeV can be accounted for by the opening up of WW decay modes. The radion 
branching ratio to a pair of gluons remains almost unchanged with mass, once the WW threshold 
is crossed. For non-zero value of the mixing parameter the branching ratio rises with the radion 
mass. This is in sharp contrast with the higgs case. Higgs branching ratio to gluons in presence of 
mixing tends to decrease from its SM value for heavy h' . 
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Figure 2: Radion and higgs branching ratios for £ = 1/6 and £ = for = 1 TeV. Whenever 
required we fix the radion and higgs mass parameters at 150 GeV. 

What should we look for at a linear collider? The lesson we have learnt from previous discussion 
is that the radion coupling to a pair of gluons is enhanced due to the trace anomaly. This has been 
exploited in previous studies on radion production at hadronic colliders. The aim of the present 
work is to see the feasibility of next generation high energy e + e~ collider to differentiate between the 
Randall-Sundrum type models with radions and SM. As radion to photon photon coupling is also 
enhanced slightly (due to the modest running of QED (3- function) , possibilities of producing radions 
in photon photon collision (with back-scattered photons, using laser and high-energy electron or 
positron beams) have been considered [13]. This cross-section is not very impressive. Moreover, 
SM higgs has same kind of production and decay channel with almost competing strength. Thus, 
this production mode may not be suitable to differentiate between higgs and radion. This is 
why we turn to the process e + e _ — > vi>(/)(—> gg). We consider both the WW fusion, which is 
the dominant production mechanism for the linear collider cm. energies, and radion strahlung, 
e+e~ — > Z* — > Zip. Two jets in the final state are easier to deal with at an e + e~ machine than 
in a hadron collider. At the same time, radion branching ratio to gluons does not fall steadily 
with its mass apart from a sudden decrease around 2W/2Z threshold. This is related to the fact 
that radion gluon gluon coupling is independent of radion mass while for higgs the corresponding 
coupling decreases rather sharply with increased higgs mass. This is crucial in our analysis in the 
sense that unlike the higgs, gluon gluon branching ratio of radion can be substantial at high radion 
mass. 

Moreover, for £ = 0, radion production via WW fusion is suppressed due to the factor (v/A^) 
in the coupling, with respect to higgs production in the same channel. Thus we have to be careful 
to choose such a decay channel for radion that this suppression factor can be overcome. This leads 



us to choose the radion decay to a pair of gluons. Even for £ / 0, we can see from fig. 1 that for 
a wide (allowed) range of £ and radion mass, the radion coupling to a pair of W's is larger than 
the same for higgs. This in turn implies the higher radion cross-section even at = 1 TeV. On 
the other hand, when 4>' — W — W coupling is small (around £ = 1/6), radion branching ratio to 
gluon gluon becomes big. This compensates the suppression in production. This already proves 
the efficacy of this particular channel. 

We already mentioned that in the limit £ — ► 0, 4> can be identified with the radion. In this 
work we will discuss the detection of <j>' through gg decay. We will not present the other decay 
branching ratios of <f>' here. These can be obtained from our earlier works [7]. 

We have calculated the cross-section for e + e~ — > vvtf (all three i/-flavours are added appro- 
priately) and multiply by 4>' branching ratio to a pair of gluons. SM higgs boson has this decay 
mode as well and produce similar kind of a signal. On the other hand, higgs decay rate to gluons 
is suppressed w.r.t radions as discussed earlier. For the purpose of comparison we have also cal- 
culated the two-jet + missing energy yield in e + e~ collision via SM higgs production and decay. 
For low mass (around 100 GeV), the SM higgs decays dominantly to a pair of b-quarks, thus also 
producing two jets. If we assume a good b-jet identification/discrimination at a linear collider and 
veto any b-jet, this gives us a handle to discriminate the SM higgs from radion (around 100 GeV 
higgs/radion mass). We have also estimated the cross-section of two jets + missing momentum 
final state coming from e + e~ — ► vv (^f*,Z*,Z — ►) qq. The strategy is to compare the invariant 
mass distribution of the jet pair. For radion one expects to have a peak around the radion mass 
we are interested in. On the other hand, the SM mass distribution of the jet pair has a continuum 
apart from a bump around Z-mass due to the on-shell Z-production. Two-jet mass distribution 
also peaks at low mass (of two jets), corresponding to the soft singularity associated with the qq- 
pairs fragmented from a soft photon. We have used the following kinematic cuts on signal and SM 
background. 

• for two jets \pr\ > 15 GeV and \rfj\ < 3.0. 

• Ar,, ( + A0g) > 0.7. 

• fa > 15 GeV. 

There is one more contrasting feature of signal and background which helps us to reduce the 
background rate without affecting the signal too much. There is a dominant part of the background 
in which the neutrinos are coming from an on-shell Z. This can be easily seen from the missing mass 
distribution of signal and background as plotted in fig. 3. Around 90 GeV, the signal distribution 
also shows a bump, corresponding to the radion production in association with a real Z. We can 
eliminate a part of the background by imposing a further cut on the missing mass: m m i s > 120 
GeV. However, the dominant part of the background remains, and it has similar topology with the 
signal events. 




Figure 3: Missing mass distribution of signal and background. 



In fig. 4 we have plotted the number of events against the two-jet invariant mass assuming an 
integrated luminosity of 500 fb _1 . For radion production the two jet invariant mass corresponds 
to the radion mass. When calculating the scatter plots, we have fixed the higgs mass parameter 
in the Lagrangian, = 150 GeV, and we have varied the radion mass parameter in the range 
50 GeV< < 800 GeV, as well as the mixing parameter in the range —0.75 < £ < 0.59. For 
the sake of comparison we have also plotted the corresponding numbers for the SM higgs boson 
(dotted lines), and for the radion in the nonmixed case, £ = (dashed lines). The upper histogram 
line shows the background coming from the processes discussed above. For the detection the ratio 
of the signal S to the background B should be S/yJ~B > 5. The lower histogram line represents 
h\f~B. Thus a 5a detection can be done if the parameter space point, marked by + is above the 
lower histogram line. As expected, the cross-sections for both signal and background grow with the 
center-of mass energy. The sudden decrease of the radion cross-section around =160 GeV is a 
reflection of decrement of <j/ — > gg branching ratio due to the opening up of WW decay channel. 
A local peak in the SM higgs contribution (only seen in \/s = 3 TeV panel), around 350 GeV of 
mass of the jet pair (same as the higgs mass) is due to the opening up of the top-pair threshold to 
the h — g — g coupling, via top loop. 

We can estimate the detection possibilities with radion VEVs other than 1 TeV, by noticing 
that the cross-sections behave approximately as 1/A^. This is evident from eqn. 6. The <j>'gg 
coupling also depends on A^, but the branching ratio of 4>' — * gg is almost insensitive to radion 
vev. Note that, apart from an overall in the couplings, quantities like a\2,0 etc. also depends 
on this vev. Thus it is easy to conclude that for A^ = 3 TeV, a very small region of the parameter 
space (mainly for < 200 GeV)would still be available with y/s = 1 TeV and a slightly larger 
region with \fs = 3 TeV. = 5 TeV or above cannot be detected by any of the studied center of 
mass energies. 



We have not separately presented the numbers for h! . However, one can see from the general 
arguments that b! production with a neutrino pair and the consequent decay to gg, is suppressed 
with respect to the 4>'. 

It is evident from fig. 1 that h! production via WW fusion is less sensitive to £, than <p' . The 
production rate is always suppressed w.r.t the (/>', except for £ values close to g. For £ values 
close to 1/6, 4>' branching ratio to gluons is orders of magnitude larger than the same for h! . We 
have explicitly checked that number of 2 jet + missing energy events coming form h' production 
and decay can only be above the 5a fluctuation of the SM background when higgs mass is below 
the WW threshold. Thus, in general one may see two resonances below the WW threshold, but 
presence of a single resonance (in gluon gluon invariant mass distribution) with mass greater than 
200 GeV, definitely points towards the presence of extra scalar like radion. 

It is evident from the figures 4 that radion-higgs mixing has a very positive effect on the radion 
search. As for example, with 800 GeV and 1 TeV center-of-mass energy, our proposed signal can 
test radion mass up to say 200 GeV in the no-mixing case. When the mixing is put on, this mass 
reach is certainly improved for favorable cases. For all the examples in fig. 4, the radion cross- 
section is well above the higgs cross-section for £ = case, while non-zero values of £ can push the 
cross-section in both ways with respect to the no-mixing case. 

The above three figures can be nicely summarised in a single plot, depicting the regions (in 
£ — rriff,' plane) in which the proposed signal is significantly above the background. In fig. 5, we 
have identified the regions where the numbers of signal events are more than five times the square 
root of the background. The negative values of mixing seem to be favorable to probe heavier 
radions for lower center of mass energies. In the shaded regions marked with a and b along with 
white region, the signal strength is higher than the 5a fluctuation of the background. One can 
claim a discovery of radions in these regions of the parameter space. Regions marked with a can 
be probed by a 800 GeV linear collider and regions marked with b by a 1 TeV linear collider. To 
probe the white regions, a 3 TeV linear collider is required. In the shaded region marked with c, 
(due to conspiracy of either the radion coupling to WW/ZZ or to gg) signal is too weak compared 
with the background even at an e + e~ center of mass energy of 3 TeV. In fact, in the shaded region 
(marked c) parallel to the m<y axis, around £ = 1/6, the radion production is too small. This 
can be accounted for by the small coupling of the radion to WW/ZZ as evident from fig. 1. A 
larger radion branching ratio to gluons in this region cannot compensate to yield a signal strength 
comparable with background. In the dark wedge shaped regions, no physical values for radion 
masses (m^) are possible for any input parameters. This can be explained from the definition of 
mixing angle 9 (see eqn.4) in terms of other input parameters. 

To summarise, we have investigated the possible signature of radions in Randall-Sundrum sce- 
nario at future e + e~ colliders. If nature chooses the Randall-Sundrum type of geometry of space 
and time, then radion, the lowest lying gravitational excitation in this scenario, can be explored at 
the hadronic collider like LHC. Being a graviscalar, it couples to the SM particles via the trace of 



the energy-momentum tensor. This implies that the couplings are very similar to that of the SM 
higgs boson. One important difference from higgs is the gluon gluon coupling, which is enhanced 
for radions due to conformal anomaly. It is important to look not only for the possible signatures 
for this scalar particle at future e + e~ machines, but also to differentiate this from the SM higgs 
boson. We propose to look for the process in which radion is dominantly produced via WW fusion 
in e + e~ collision along with neutrinos and subsequently decays to a pair of gluons. This in the final 
state produces 2 jet + missing momentum signature, with a peak in two jet mass distribution. We 
have also estimated the contribution to this final state from SM higgs boson and also from other 
SM backgrounds. Though at an e + e~ collider with yfs = 500 GeV the detection is not obvious, at 
y/s = 800 GeV and above one can expect 5a effect in exploring the radion and differentiating it 
from the SM higgs boson. 
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Figure 4: Number of signal (+ for ( / and dashed line is for £ = 0), SM higgs (dotted line) 
and background (solid histogram) events as a function of the invariant mass of the jets. The 
upper histogram hows the actual number of background events. The lower histogram shows the 5er 
fluctuation of background. Points (representing the signal from higgs and radion) above the lower 
histogram can be explored at 5a. Whenever required we fix the higgs mass parameter at 150 GeV. 
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Figure 5: Radion mass reach as a function of the mixing parameter £. Radion masses, corresponding 
to the dark wedge-shaped regions on both sides of £ = 0, are not allowed. Regions (triangular area 
close to £ = and in the left side) marked with a, can be probed at 5a level with e + e~ center of 
mass energy 800 GeV. Region marked with b along with a can be probed at 5a level with e + e~ 
center of mass energy 1 TeV. The white regions along with a and b can be probed with 3 TeV e + e~ 
center of mass energy. In the shaded regions marked with c (note that there is a small region c 
close to £ = but for £ < 0) either the radion production is too small or the gg decay branching 
ratio is too small. These regions cannot be probed. 
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